Dual formulations of non-abelian spin models: 
local representation and low-temperature 

asymptotics 

O. Borisenko^, V. Kushnir^ 

N.N.Bogolyubov Institute for Theoretical Physics, National Academy 
of Sciences of Ukraine, 03143 Kiev, Ukraine 

Abstract 

Non-abelian lattice spin models with symmetry group SU{N) or U (N) 
can be formulated in terms of link variables which are subject to the Bianchi 
constraints. Using this representation we derive exact and local dual formu- 
lation for the partition function of such models on a cubic lattice in arbitrary 
dimension D. Locality means that the dual action is given by a sum over 
some subset of hypercubes of the dual lattice and the interaction between 
dual variables ranges over one given hypercube. Dual variables are in general 
discrete-valued and live on {D — 2)-cell of the dual lattice, in close anal- 
ogy with the XY model. We use our construction to study in details the 
dual of SU{2) principal chiral model in two dimensions. We give dual ex- 
pressions also for two-point correlation function in arbitrary representation 
and for the free energy of defects. Leading terms of the asymptotic expan- 
sion of the dual Boltzmann factor are computed and it is proven that at low 
temperatures it converges to a certain Gaussian distribution uniformly in all 
fluctuations of dual variables. This result enables us to define the semiclassical 
limit of the dual formulation and to determine an analog of the vortex-spin- 
wave representation for the partition function. Such representation is used 
to extract leading perturbative contribution to the correlation function which 
shows power-like decay at weak coupling. We also present some analytical ev- 
idences that the low-temperature limit of the dual formulation is completely 
described by /5'0(2)-like approximation of SU{2) matrix elements. 



1 Introduction 

Two-dimensional {2D) lattice spin models, besides being interesting on their own 
right as models of ferromagnets can be viewed as mathematically well defined 
scheme of the nonperturbative regularization of quantum continuum field theories 
like SU{N) principal chiral model or 0{N) nonlinear sigma model. Being a dis- 
crete theory lattice models can be simulated on computers by Monte-Carlo method, 
and this is nowadays the most important tool of obtaining physical results near the 
continuum limit. Among major analytical methods one could mention the strong 
coupling expansion and the perturbation theory (PT). The PT is essentially the 
only analytical tool which provides systematic expansion of different physical quan- 
tities in the region of weak coupling, i.e. in the region relevant for the construction 
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of the continuum limit. However, at the best PT is only applicable for studying 
short- distance quantities and the most important and interesting phenomena, like 
the mass gap generation cannot be described in its frameworks. When attempting to 
go beyond PT in the region of weak bare coupling, one runs into various mathemat- 
ical problems, the most important being an absence of any analytical control over 
exponentially small contributions and, hence the absence of any reliable method to 
study long-distance physics where nonperturbative effects dominate even if they are 
exponentially suppressed. These facts impelled people to look for different, though 
equivalent representations which would allow to study the long-distance physics. 
One of such popular and promising representations is known as a dual formulation 
and is based on a certain non-classical change of variables. This formulation deals 
with dual lattice and appears to be very fruitfuU for abelian models PP. Namely, 
within such formulation 1) have been obtained various important analytical results 
on the long-distance dynamics of abelian spin fields and 2) the topological structure 
of the vacuum is transparent and can be more easily studied. 

In the context relevant to this paper we would like to mention the dual of 
the abelian XY model pQ which has been used to prove the existence of a soft 
phase at low temperatures with power-like decay of the correlation function in two- 
dimensional XY model |2j. In this case the dual of the XY model is a local theory 
for certain discrete variables. The conventional dual transformations |SI for non- 
abelian models are perfectly well defined in mathematical terms and also lead to a 
local dual theory for integers which label irreducible representations of the symme- 
try group. Nevertheless, these transformations are not complete if one compares to 
the abelian case. First of all, the resulting dual variables are not independent but 
are subject to some constraints known as triangular conditions, and as such they 
cannot really be associated with elements of a dual lattice (of course, introduction 
of the dual lattice is not the crucial point but rather matter of convenience; e.g., 
in the abelian case one can work with dual variables on the original lattice). It 
is to say, however that it is not really clear what are true independent degrees of 
freedom of the duals of non-abelian models. Secondly, such formulation is so math- 
ematically involved that today it is even far from obvious if it can be useful for 
any kind of the study of the model, except for the strong coupling expansion. In 
particular, the positivity of the Boltzmann weight remains unclear (see discussion in 
[2). Even the precise definition of what can be termed the dual Boltzmann weight 
is not really obvious, so the possibility of numerical simulations is at present rather 
uncertain. Most importantly, however is that it is not clear how one could proceed 
in an analytical study of the model at low-temperatures^. 

On the other hand, there exists representation of SU{N) and U{N) spin models 
in terms of link variables and this representation can be formulated directly on 
the dual lattice. It is a first goal of the present paper to use the link representation 
for the derivation of an exact dual formulation of non-abelian spin models. Here we 
give two such formulations which appear to be quite different from the conventional 
formulation mentioned above. In our opinion the most essential advantages of our 

■^This is, of course personal opinion of the authors. After some attempts to understand low- 
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formulations are that 1) one of them can be definitely used for the Monte-Carlo 
simulations in 2D and 2) they are much more suitable for an analytical investiga- 
tion of the model in the low-temperature region. The second fact follows from the 
properties of the link formulation and we refer for the details of why it is so to our 
papers [H]. In the last of those papers we have already presented a dual of 2D 
SU (2) spin model and proposed approximate representation for the dual partition 
function at low temperatures. Another important feature of one of our formulations 
is that it carries a close analogy with the abelian U{1) model. For instance, for 
2D SU (2) model there are 3 independent degrees of freedom per (dual) site. At 
low temperatures it possesses an analog of the vortex-spin-wave representation, etc. 
The deficiency of our approach is that it can be straightforwardly applied only to 
models in which spins are elements of the Lie group. We do not know any obvious 
generalization of link formulation to, say, 0{N) nonlinear sigma models. 

Low-temperature properties of 2D non-abelian models are crucial for construc- 
tion of their continuum limit. It is commonly recognized that models possess no 
phase transition, correlation function has exponential decay at any coupling and 
models are asymptotically free. Despite being more than twenty five years old this 
expectation has not been proven rigorously. On the contrary, certain percolation 
theory arguments suggest that all non-abelian models have soft low-temperature 
phase with power-like decay of the correlation function j7j. It is thus another im- 
portant motivation of the present investigation, namely to get deeper insight into 
the nature of the mass gap in non-abelian spin models. For example, in many pa- 
pers devoted to 2D non-abelian models it is written that "there is a nonperturbative 
mass gap generation at arbitrarily small couplings". It is not really clear, however 
what is precise meaning of this "nonperturbative generation" . It cannot be a sim- 
ple consequence of the link decorrelation which happens in ID models. Then, one 
could ask if this "nonperturbative generation" follows from the existence of some 
non-trivial background of defects like vortices of the XY model or is due to the 
strong but smooth disorder of non-abelian spins, e.g. like center vortices |H]. As 
is well known, the dual formulation of abelian models have been extremely useful 
in clarifying all these important physical problems |21 El Cn]- We think that the 
dual formulation given here is a good starting point for obtaining reliable analytical 
results in the low-temperature limit of non-abelian models. This is our next goal 
to develop a technique within dual formulation of non-abelian models which would 
allow to investigate them in the limit of the weak coupling. In this context we study 
here the following two approaches on the simplest example of SU{2) spin model 

• We derive an asymptotic expansion of the dual Boltzmann weight when /3 — > cxo 
uniformly valid in fluctuations of the dual variables. The same procedure is 
also done for the two-point correlation function. As will be seen, at low tem- 
peratures the Boltzmann weight converges to a certain Gaussian ensemble - 
the fact which ensures a possibility to define an analog of the vortex-spin- wave 
representation for the partition and correlation functions in the semiclassical 
limit. This step provides a major simplification of the model as /? — cx) but 
even in this case it cannot be yet solved exactly. The formulation obtained 
permits very simple calculation of the leading perturbative contribution to the 



correlation function. This contribution results in the power-like decay of the 
correlation function. 



• The second approach consists in replacing SU{2) matrix elements by ISO{2) 
ones in the vicinity of the identity element of SU{2). We shall give a proof 
that such replacement is valid at low temperatures, compute the asymptotics 
of the dual Boltzmann weight and present simple mimcrical evidence that prac- 
tically for all configurations this way of calculations produces very reasonable 
approximation for the original Boltzmann weight. 

This paper is organized as follows. In the section 2 we review briefly the conven- 
tional dual transformations. After introduction of the link formulation we develop 
original approach to the dual transformations and present two forms of the duals 
of the principal chiral model. In the section 3 we investigate the low-temperature 
properties of the dual Boltzmann weights of SU{2) model in two dimensions along 
the lines described above. Some technical details of this investigation are given in 
the appendix. In the section 4 we summarize our results and outline some per- 
spectives for future investigations. Also, we calculate here the leading perturbative 
contribution to the correlation function. 

2 Principal chiral models and their duals 

We work on a /^-dimensional hypercubic lattice A E Z'^ with lattice spacing a = 1 
and a linear extension L. We impose either free or periodic boundary conditions 
(BC). Let G = U{N), SU{N); U{x) e G and DU{x) denote the Haar measure on G. 
Xr{U) and d{r) will denote the character and dimension of irreducible representation 
{r} of G, correspondingly. We treat models with local interactions, i.e. interactions 
between nearest neighbours. Let H[U] be a real invariant function on G such that 



exist. The partition function (PF) of the principal chiral model with the symmetry 
group G X G is defined as 



H[U] I < H{I) 



(1) 



for all U and coefficients of the character expansion of exp{PH) 





MP)- llDU{x)exp J]i/ [t/(x)C/t(x + e„)] 



(3) 




2.1 Standard dual transformations 



Conventional dual transformations can be defined as a sequence of transformations 
consisting of the following steps: 



1. Fourier expansion of the Boltzmann weight exp(if). This is an essentially 
character expansion on the group. 

2. Exact integration over original degrees of freedom. As a result one obtains a 
set of constraints on the summation variables (which label representations of 
the group G and matrix elements of U{x)) in the character expansion. 

3. Solution of the constraints in terms of new dual variables. 



Performing the first step one finds for the PF 
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where / = (x, e^) denotes links of the lattice and Q{x) is a group integral. Let 
D^'^[f/] be a matrix element of U in representation r. Then Q{x) can be written as 
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The second step in the most general settings (for arbitrary graph and for all G) has 
been accomplished in jS]. Through the rest of this subsection we consider only the 
simplest cases of G = U{1) and G = SU{2) in two dimensions. These examples 
already show a marked difference between abelian and non-abelian models, on the 
one hand. On the other hand, we want to have relatively simple formulae which can 
be compared with our formulations. For the U{1) model we have 



Q{x) = 6 [r{x, en) - r{x - €„, €„)] 



(6) 



>n=l 



where 6{x) is the Kronecker delta. In SU{N) case the result for Q{x) in j^] is 
expressed in terms of the Haar intertwiners. We re-state this result in more familiar 
form expanding Q{x) in the Clebsch-Gordan (CG) series 
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where ri = r(a;,ei), r2 = r(a;,e2), = r{x — ei,ei), = r{x — 62,62) and similar 
notations are used for mi and fcj. 

The third step can be readily accomplished for the U{1) model. We skip all the 
details of the derivation which are well known and can be found, for example in 
The result for the PF reads 
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where the product runs over all links of the dual lattice A and the dual Boltzmann 
weight is given by 

^ S^^(r(x) - r{x + e„); P) = ^ exp + i{t{x) - t{x + e„))] , 

Jo ^TT 

(9) 

which in the case of the standard choice H = cose/) leads to 

^0^ (l) = Ir{x)-r{x+e^)iP) (10) 

with - the modified Bessel function. 

Unfortunately, situation is much more difficult for non-abelian case. Q{x) for 
SU (2) contains not just one but 4 conditions: 2 Kronecker delta's which constrain 
magnetic numbers and ki 

mi + 1712 = + = s , ^3 + 1714 = ki + k2 = t (H) 

and 2 triangular conditions on representations 

\ ri - r2 \< J < ri + r2 , \ rs - \< J < r-^ + . (12) 

We are not aware of any attempts in the literature to resolve these constraints 
simultaneously for all lattice sites. Of course, this does not cause much problems 
for the strong coupling expansion where only lowest values of r{l) contribute. But 
situation becomes hopeless in the low temperature region where probably all r{l) 
become relevant. Eventually, this reduces to the problem of finding the appropriate 
summation technique when (3 —>■ oo which we could not solve. Moreover, Q{x) as 
defined in ((7j) is not strictly positive quantity: on a number of configurations it is 
negative. We do not know how to perform resummation over the magnetic numbers 
to get positive Boltzmann weight. Most probably, only complete summation over all 
magnetic numbers will produce positive quantity. One of the possible ways of how 
to proceed further is to go over to the triangular lattices identifying representations 
with the length of the sides of triangles and magnetic numbers with the orientation of 
the triangles, and we have undertaken some attempts to go along this line. Detailed 
description of this procedure is however beyond the scope of the present paper. 



2.2 Link formulation of principal chiral models 

We turn now to a different approach to duality transformations based on so-called 
link representation for partition and correlation functions. That is, the three steps 
formulated in the beginning of the previous subsection are replaced by the following 
ones 

1. Change of variables in the partition function Q 

V{1) = U{x)U\x + en) . (13) 



Integration over original degrees of freedom. This integration generates a set 
of constraints on the link matrices V(l) known as the Bianchi identities. 



2. Implementation of constraints into the partition function by making use of 
the invariant delta-function on the group. This step introduces new degrees 
of freedom which can be associated with plaquettes of the original lattice and 
which label the irreducible representations of the group G. At this step one can 
go over to the dual lattice so that the new variables belong to the elements of 
the dual lattice. Integration over link variables. These two steps give already 
some dual version of the principal chiral model. 

3. Before integration over link variables one can decouple summations over ma- 
trix indices employing certain orthogonality relations for the CG coefficients 
(or for some other objects) of the group and multiplying the group matrices 
entering Bianchi identity in some specially ordered way. Though details de- 
pend on the dimension of the lattice the result is that only independent dual 
variables remain in the theory and their interaction is confined to a subset of 
all hypercubes of the dual lattice. 

The first step of this program has been accomplished twenty years ago in Q . We 
therefore restrict ourselves only to the description of the result. The PF can be 
exactly reformulated in terms of link variables as (with periodic BC) 



/3 5^i/[V(/)] 



with Jacobian J[V] given by 

J[V] = n J[V{p)] n AH{n)] 

peA n=l 

where J[V"(p)] and J[H{n)] are given by 
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J[V{P)] = E^W^^'(n^(^)) ' 
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J[H{n)] = 5^rf(r)x. ( n^(^«)) ' ^ = (^.e„). 
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(16) 



(17) 



Yip is a product over all plaquettes of lattice A, the sum over {r} is sum over 
all representations of G = U{N), SU{N), d{r) = Xr{I) is the dimension of r-th 
representation. The group character Xr depends on a product of the link matrices 
^(0 = yn{x) along a closed path (plaquette in our case): 

W V{1) = Vn{x)V^{x + en)V},{x + e^)Vi{x) . (18) 

The expression (irXrdliep ^(0) is the invariant group delta-function which re- 
flects the fact that the product of U{x)U'^{x + e„) around plaquette equals I 



\{V{l) = I. 



(19) 



This constraint is called local Bianchi identity while the constraint on holonomy 
operators 



: / (20) 

ln = l 

is usually called global Bianchi identity. The only solution of two constraints p9|) 
and ipUl) is a pure gauge (fT^ which restores the equivalence of the standard and 
the link formulations. On the lattice with free BC the global constraint should be 
omitted. Furthermore, in the thermodynamic limit (TL) the contribution from the 
constraints on the holonomy operators vanishes very fast. In the perturbation theory 
in 2D it vanishes like inverse powers of L [S], and similar behaviour is expected 
for D > 2. Since we are eventually interested in the TL we neglect these global 
constraints (in fact, the presence of global constraints does not cause any principal 
problems as will be seen from the procedure described below but their absence makes 
the technical details somewhat simpler). Thus, for both free and periodic BC we 
work with the PF 
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The two-point correlation function in representation j takes the form 

r,{x,y) = (d-\j)xj I n ^(0) ) > (22) 

\ \leCxy J I 

where Cxy is some path connecting points x and y and W{1) = V{1) if along the path 
the link / points positive direction and W{1) = V^^(/), otherwise. For more details on 
the link formulation we refer the reader to our paper where we have developed 
a weak coupling expansion for SU {N) spin models using the link representation. 



2.3 Dual of the link formulation: I 



Here we are going to accomplish the step 2 described in the previous subsection. The 
idea that the true dual variables in the link formulation are the Fourier conjugate to 
the local Bianchi identity has been proposed also in |lj but the formulae have been 
explicitely written only for abelian U{1) model in 2D. We thus want to reformulate 
the model !^21\i on the dual lattice only in terms of discrete variables that are in our 
case representations {r{p)} and magnetic quantum numbers mi{p). We shall use the 
following conventions for the dual lattice, fl^ will denote an object dual to site x, 
i.e. D-cell of the dual lattice, Qi - {D — l)-cell of the dual lattice, Qp - {D — 2)-cell 
of the dual lattice and so on. By definition, object dual to k-th cell is {D — A;)-cell. 
As follows from (PT|) and from the definition of the character which enters the group 
delta-function in (fTB|) XriV) = ^„ V^™ on the dual lattice the PF may be written 
as 
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(23) 



Due to the trace there are 4 variables rrii at each Qp, thus i = 1,2,3,4. The dual 
weight EqIQi) is given by the following one-link integral 



D-l 



(24) 



where V^"" is a matrix element of r-th representation. {D — 2)-cells Qp and Qp share 
{D — l)-cell fli. Ylu runs over {D — 1) such pairs. 

Similar form for the two-point correlation function in the representation j reads 
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where sn+i = si and the link integral on Qi e Cxy is 
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(26) 



Here, C^^y is a path dual to the path Cxy between points x, y on the original lattice. 

To vizualize last formulae let us give expression for the PF in two dimensions 
with G = SU{2). In 2D we have = P, = ^ and Qp = x. Placing as usually 
dual sites in the centers of original plaquettes we obtain 
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(27) 



where the dual weight 

So(0 = Eq {r{x),mi{x),mi+i{x);r{x + en), ti{x + en), U+i{x + en); P) (28) 
becomes 

-^(^l'^ _ J (IYqI^HIV] yruiix) nii+iix) y]^ ti(a;+e„) ti+i(x+e„) ^29') 

It is very easy to make an integration in expanding the result into the CG 
series. One then finds the following representation for the dual weight of SU (2) 
model ^ 

^ 2r I 1 '^'"^^'^^ '^Ifm^ Jk Clfm2 Jk ' (30) 
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where we have denoted ri = r{x), r2 = r{x + en), mi = mi{x) and so on. For the 
standard choice H[V] = Xi/2(V") the coefficients C[J] can be computed explicitely 



C[J] 
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For the correlation function in Eq. (j25p using the CG expansion one gets 

Jaia2 

Last expressions can be compared with the dual of the U{1) model given in (jH)) and 
with the standard dual formulation of SU{2) model described by formulae ^ and 
d?!). Certainly, the dual weight for SU{2) model Ho(/) as given in (j^ - lj^ is very 
close to what we have for the dual weight in U{1) model: 1) dual variables, i.e. 
representations r{x) and magnetic numbers mi{x), live on sites of the dual lattice; 
2) the full Gibbs measure is factorized into product of dual weights So(/) over links 
of the dual lattice. Unlike the standard dual formulation there is no complicated 
triangular conditions on representations: triangular constraint on J in ()30p is only 
one-link problem which does not lead to any significant complications. Nevertheless, 
there is also a difference from the abelian case. As follows from the properties of 
the CG coefficients we have one constraint on the magnetic numbers on every link, 
namely summation over k in (jHUj) produces the following condition 

mi{x) - mi+i{x) = ti+i{x + en) - ti{x + en) . (33) 

This shows that not all dual variables are in fact independent, and this is the main 
difference from the U{1) model. There is no simple way to make a change of variables 
such that one gets only independent variables and keeps the locality of interaction. 
In the next subsection we present a different way of the construction of the dual 
theory from the link formulation in such a way that the dual weight depends only 
on independent dual variables. 

We finish this subsection with the brief description of some important features 
of So(0: 

• As follows from the properties of the coefficients of the expansion C[J] the 
series in J in ()3U|) gives directly the strong coupling expansion of the model 
written in closed and compact form. Much less trivial task is to get weak 
coupling expansion for Ho(/) since all J in the series ()30|) become relevant. 

• On all configurations {rj,mj,nj} So(/) is strictly positive So(/) > 0. Though 
we could not prove it rigorously this claim is supported by the following facts: 
1) the first term in the strong coupling expansion is strictly positive, thus at 
sufficiently small j3 where the series converges very fast So(/) is positive; 2) 
the leading term of the asymptotic expansion of So(/) at large /? is strictly 
positive on all configurations; 3) numerical computations of Ho(/) on a number 
of configurations and in a wide region of /3 also support this conclusion (we have 
checked this statement using Mathematica on more than 10^ configurations). 
If So(/) > on all configurations this gives a chance for a numerical Monte- 
Carlo simulations of the dual model. 

• The dominant contribution to Ho(/) at large /3 comes from diagonal compo- 
nents of the rotation matrices, non-diagonal contribution is suppressed roughly 
as Km — n)\]~^. This is, of course a conseauence of the fact that when B ^ oo 



the link matrix performs only small fluctuations around unity. In turn, this 
property gives a possibility to compute low-temperature asymptotic expansion 
of Ho(/). This will be subject of the next section. 



2.4 Dual of the link formulation: II 

The main idea of the following approach is to introduce the independent dual vari- 
ables from the beginning, namely to find a representation for the group invariant 
delta-function in terms of quantities which would be free of any conditions. We 
explain the realization of this idea for 2D SU (2) model and then show how it can 
be extended to other groups and to higher dimensions. 

We start with the representation for SU{2) matrix elements in which the depen- 
dence on the magnetic numbers enters only through the Clebsch-Gordan coefficients 

/^aa 4 

Vri^, 0,4>) = Y1 ^rl xk Ur'^iu;, 9, 4>) , (34) 
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where {u, 6, (p) is expressed through the spherical harmonics Yxk{0, 0) and gen- 
eralized characters Xa(^) '^^ SU{2) 



U^\u,e,cP) = (-O'^^y^f^ Y,,(^,0) XaM • (35) 

Generalized character of rank A in representation r can be defined through the 
relation 

r 

XU^)=^' E e"™"^™Ao- (36) 
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We remind that the invariant measure in this parameterization is 



4vr2 7o 2 
and the fundamental trace becomes 
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Tr\/i/2 = 2cos| . (38) 
The basic formula which we need is the following representation for SU (2) characters 



(2r + 1) r r 

XriyiV2) = ' d^ dasma^{r,^,a-V{)^{r,^,a-V2) , (39) 

$(r, yp, a- V) = ^A(t) X'x H , (40) 

A=0 ^ 

^In what follows we use notations and conventions of Ref._n^ for the Clebsch-Gordan coefh- 
cients, spherical harmonics, etc. 



where Px{t) is the Legendre polynomial and 



t = cos 9 cos a + sin 6* sin a cos(0 — ip) . (41) 

This representation can be easily proved if one uses orthogonality relations for the 
spherical harmonics and then the addition theorem for the generalized characters. 

As the next step, we divide all sites of the dual lattice into sets of even and odd 
sites. Site x = {xi,X2) is even if both xi and X2 are even or odd, and is odd if one 
of Xi is odd. Thus, all even sites are surrounded by odd sites and vice-versa. Using 
cyclic properties of the trace we can couple the four matrices entering local Bianchi 
identities in two different ways, namely 

Xr{V,) = XrlVuVs,] , (42) 

if the site x is even and 

XrW = Xr[V23V,^] , (43) 

if the site x is odd. We introduced here obvious notations 
Vi2 = V{h)V{l2), Vs, = V\h)V\U), V23 = V{l2)V\k), = V\U)V{k). (44) 




Figure 1: Dual even plaquette and coupling of matrices in group delta's. = V{li), 
sites X and x + ei + 62 are even, while x + ei and x + 62 are odd. 

With help of (jnni) and (jUl), (0 the Jacobian J[V] defined in and ^ can 
be written on the dual lattice as 

j[v] = n j[yi^)] n ^[^(^)] ' (45) 

^.even ^odd 

where in even sites 

■Ji^x) = 22 A — 1 '^^ I sin a $(r, 1/12) $(r, 1/34) , (46) 

r=0,l/2,.. ^ 

while for odd sites 



^ 7 / / c?asina$(r,<^,a;F23) $(r, a; V^4i) • (47) 



». — n 1 /n 



A set of variables ^(x) = {r{x),a{x),(f{x)} is chosen to be dual variables at site 
X. Finally, we decompose all plaquettes of the dual lattice into sets of even and 
odd plaquettes in such a way that all even plaquettes are surrounded by odd ones 
and vice-versa. Due to the special coupling of matrices as described above, the 
integration over link matrices is now coupled only within every, for example even 
plaquette. This is shown in FigHJ 

Substituting (j46p and (j47j) in (jl4p . re-denoting notations for links I G p uniformly 
for all even plaquettes as shown in FigQwe obtain dual expression for the partition 
function 



^SU{2) 
A 



n 



x-gA 



(2r(x) + 1 

r{x)=0,l/2,... 



'o 



da{x) sina(x) 







X JJ 5o[e(x),e(x + ei),e(x + e2),e(a; + ei + e2); /3] 
p,even 



(45 



where the dual Boltzmann weight 

5o(p) = 5o [e(x), e(x + ei), e(x + 62), e(a: + 61 + 62); /3] (49) 

is given by 

B,{p) = /"n^^*exp|/?5^/7[Vj| $(e(x);ViV2) (50) 

1=1 L 1=1 J 

X $(e(x + 62); V2V3) ^{iix + ei + e^y.vM) + ei); V,Vl) . 

Dual expression for the correlation function can be easily re-constructed from Eq.(|22j) 
and last formulae. Consider, for simplicity the shortest path between points x = 
(0,0) and y = (0,-R). The corresponding dual path is a set of links which point 
direction ci. Let \ x — y \= Rhe even. Since every link belongs to one and only one 
even plaquette we obtain 



2j + l 11 



zeA 



(2r(x) + l)^y ' J da{x)sina{x) 



3 



r-(x)=0,l/2 
3 



X 



P' si=-3 SR=-j p(C) 



(51) 



where p' denotes a set of even plaquettes which do not contain links I G C^y and 
p{C) - a set of even plaquettes which contain such links. In notations of FigQ we 
have 

i?,(p(C);s,_iS,s,+i) = I'lldViexphY.HivA V;'-''\h) V;'''+' (k) (52) 

1=1 I i=i J 



x^(^(x):VM ^(((x + eo):V^V-,) + d + eo):V.^V}) ^>(^(x + e,):V^V:) . 



4^ 



Finally, the disorder operator in the dual representation can be readily obtained 
from its expression in the link formulation. Let H[—V{1)] = —H[V{1)]. Consider 
the following disorder operator 



D{x,y) = exp 



2(3J2H[-V{1)] 



leT 



(53) 



where T = {/} = [(x, ei), {x + 62, ei), {y — 62, ei), {y, ei)]. Expectation value of 
this operator can be written as a ratio of two partition functions. With help of 
the change of variables V{1) — > —V{1) this ratio can be presented as the following 
expectation value in the dual formulation 

{D{x,y)) = ((-i)^^')-^.)^ _ ^54) 

Generalization of these representations to higher dimensions is straightforward. 
In this case we couple link matrices in such a way as to decouple integration over 
them into a set of integrations within each even D-cell of the dual lattice. Then the 
partition function can be written as 



zf(^)(/?) = n 



{2r{n,) + ir r r sin a{n,)da{n,) 

, fi -1/0 •-'0 Jo 



r(f7p)=0,l/2,... 



X 



n B{hD) , (55) 
/id, even 



where the dual Boltzmann weight reads 

B{hD) = I n [dV{ni)exp{f3H[Vm}] J] m%y,l{^p)) ■ (56) 

Here, ^{flp) = {r{Qp), a{Qp), ip{Qp)} are dual variables which belong to (Z) — 2)-cells 
of the dual lattice and 7(f2p) = {u{Qp),9{Qp),^{flp)} are combined angles arising 
from multiplication of {D — l)-cell matrices (original link matrices ) which belong 
to the same plaquette trace. 

To our opinion, the dual representation presented above bears the closest analogy 
to the abelian case. For example, in two-dimensional model there are three inde- 
pendent dual variables per site as could be expected on general grounds. Further, 
as follows from the Eq. (j5U|) the interaction between dual variables is local though 
it is not exactly the same as in abelian models: in addition to nearest-neighbour 
interaction, there is an interaction over diagonals of even plaquettes. The deficiency 
of this dual representation is that the dual weight is in general complex valued. This 
makes the numerical MC computations impossible. On the other hand, we believe 
it is well suited for an analytical investigation of the model at low temperatures. 

In the described above scheme the realization of ideas leading to the local dual 
formulation relies on the special representation for SU{2) matrices ()34|) . (jH^j) . This 
representation appears to be very efficient at large (3 as we shall show in the next 



section. We do not know if similar parameterization exists for any group G. There- 
fore, for arbitrary group G we proceed as follows. Instead of (jSHI) we use unitarity 
relations for the CG coefficients of U{N) or SU{N) group |12j 



{jk} 



^Tim-i r-zm-z ^riuj rznz ~ ^mi,ni^m2,n2 (57) 



which allows to present group characters as 

Xr {V,V2) = J2 ^(^' ^' ^i) ^(^' ^' ^2) , (58) 

{Jk} 

where 

G{r, J, k;V) = ^ y^i"^^ . (59) 

mi ,m2 

We then proceed exactly as for SU{2). For the sake of simplicity let us consider the 
two-dimensional theory. Identifying ^{x) = {r{x), J{x), k{x)} as new dual variables 
the PF can be written as 

z^:^i^\f3) = ^n^[^(^)] n B{p), (60) 

where the dual Boltzmann weight takes the form 

B{p) = /"nrfV^iexp|/5^i/[V^,]| G(e(x);V^ir2) (61) 
i=i I i=i J 

X G(e(x + 62); V^Vs) Gi^x + ei + e2y,VM) ^(^(x + ei); V,Vl) . 

Substituting into the last equation one can express B{p) in terms of one-link 
integrals defined in (j^ . Since the CG coefficients of SU{N) group can be chosen 
real ^2]; the dual weight in this representation is also real. We have not studied if 
it is positive. 



3 2D SU{2): Dual weights at low temperatures 

One of the most important application of the dual formulation is the investigation of 
the low-temperature region of two-dimensional non-abelian models. The first step in 
this investigation is to establish an asymptotic expansion for dual weights at large /3. 
Clearly, it is necessary to get asymptotics uniformly valid in all fiuctuations of dual 
variables. It turns out that such asymptotics can indeed be constructed, and this 
is, in our opinion one of the most important advantages of our dual formulation. 
Calculation of the asymptotic expansion essentially relies on the fact that when 
(3 —>■ 00 the link matrix performs only small fluctuations around unity both in the 
finite volume and, most importantly in the thermodynamic limit. Note, that this is 
not the case for the original {U{x)) degrees of freedom: in the large volume limit 
their fluctuations are not bounded. 



3.1 Effective theory at low temperatures 



In this subsection we construct an effective low-temperature theory for the dual 
formulation PH|) . The corresponding dual weight is given by Eq. ()50|) . We modify it 
by introducing sources jk{l) which can be related to the correlation function 

Bj{p) = /"nrfV^exp|^[/3/7[V,]+jfe(/K,(/)]| ^a^y,V,V2) (62) 
1=1 1 1=1 J 

X + 62); ^^2^3) + ei + 62); V^) + ei); V,V^) , 



where 



cui (/) = sin ^ sin ^ cos $ , ct;2(0 = sin — sin ^ sin $ , c<j3(/) = sin ^ cos^ . 



It is well known (see, e.g. [TB]) that when /3 — >■ oo the link matrix performs only 
small fluctuations around unit matrix. Our next calculations relay on this fact. 
Let H[V] has a unique maximum at V = I. We also require that H[V] possesses 
character expansion. Then, in the vicinity V = I and in the parameterization (j34j) . 
()35|) one can always construct the following expansion 

H[Vi] = H[I] - 1 sin' ^ + 0{uj^) . (63) 

For the standard choice H\V] = TtVi/2 = 2cos^ one finds H[I] = 2 and 7 = 2. 
Next, we need uniform asymptotic expansion of the function $(^;V). We derive 
such an expansion in the appendix A. The leading term of this expansion takes the 
form 

ln<l>(^;l^) = -i sin '^WR'^ - 1 + ^(sin^ |) , (64) 

where the classical angular momentum R is introduced in (jHH). As follows from the 
semiclassical expansion pOlj) . the remainder is actually bounded as 0{l3~^) for all 
R. Therefore, corrections 0{sm' ^) and higher are in fact negligible in the limit 
/3 ^ 00 so we treat them perturbatively. Lastly, we substitute all these expressions 
into ()62|) . express combined angles u, 6, $ in terms of link angles and perform 
integration. We omit all these simple but rather cumbersome calculations and give 
only the final result. Introducing notations 

C{(3) = exp [2(3L'H[I] - 3LHn{jf3)] , (65) 



Rk{x) = ^= tkix)^/R^ix)-l , (66) 



ti{x) = sma{x) cosip{x) , t2{x) =sina;(x) sm(p{x) , ti{x) = cosa;(x) 

and making the change of variables R{x) = 2r{x) + 1 so that R{x) takes integer 
values, we present the result for the dual partition function at large f3 in the form 



y K [x) j — / aa(x) sma[x 



R(x)=l 



exp[Seff] . (67) 



An expression for the effective action Sejf is too large to be given liere in full, 
therefore we present only its expansion. Rescahng sources as jk{l) = -^jk{i) the 
effective action can be expanded as 



S, 



eff = So + ^S, + ^S, + {(3-'/') 



1/3 



(68) 



We have computed first three terms in the above expansion. Here we give 5*0 and 
^1 1 

5*0 = --^{Rkix) - Rk{x + Cn) + ijk{l)) , 



S, = e'"^' J]l^,^(x)/2fc(x) , 

X 

Wim{x) = Pi{li)pmil2) +Piih)Pmik) X is even 



(69) 
(70) 
(71) 
(72) 



Wimix) = Pi{li)pm{k) +Pi{h)Pm{l2) , if x is odd . 
We have introduced here the following notations 

Pk{li) = Rk{x) - Rk{x + ei) + iik{x,ei),pk{h) = Rk{x) - Rk{x + €2) +ijfc(x, 62), 

Pk{h) = Rk{x-ei)-Rk{x)+ijk{x-ei,ei),pk{l2) = Rk{x-e2)-Rk{x)+iik{x-e2, €2) . 

Due to the term R^{x) in the measure, summation over R{x) can be extended to 
include point R{x) = 0. We use Poisson summation formula, perfrom change of 
variables R{x) ^/^R{x) and make the following expansion 



Rkix) 



tk{x)^ll3R\x) 



tk{x)R{x) + 0{(3-') . 



(73) 



This expansion is valid when '-y(3R^{x) > 1. We call this expansion semiclassical 
approximation. It allows to define an obvious analog of the vortex-spin-wave rep- 
resentation for the PF 



C{f3)l[ 



k(X 



k=l 



(74) 



X exp 



m{x)=—oo ' 

S,ff{Rk{x)) + 2t'K^J2m{x) I ^2^1 



1/2 



X 



3.2 /5'0(2)-like approximation 

Here we investigate the dual model given by Eq.(|27p with the corresponding dual 
weight (j29|) . Our approach is to replace the matrix elements in the integrand of 
by the /5'0(2)-like matrix elements. Justification of such replacement will be 



given in the course of calculations. For SU{2) matrix elements we use the Wigner 
D-function parametrized by Euler angles 

Dl^" {a, u, 7) = e-™'^-*"^rf;""(cj) . (75) 

For the sake of simplicity we consider here only the standard action H[V] = J2n ^1%^ 
i.e. the fundamental character which in this parameterization reads 

5^Dr;2 = 2cos|cos^(a + 7) ■ (76) 
The invariant measure on the group takes the form 



dV = ^ / ^7 / da dusinuj . (77) 
lOTT Jo Jo Jo 

Substituting last expressions into ()29|) one can exactly integrate over a and 7 angles. 
This gives 

2o(0 = C™ f^c^sinu;cosa; /2(„,-t,)(2/5 cosu;) rf™^'"^(2a;)<f (2a;) . (78) 
Jo 

To get the asymptotics when /? ^ 00 we use the fact (described in the previous 
subsection) that the link matrix in this region performs small fluctuations around 
unit matrix, namely 

00 ~ 0{(3-^/^) . (79) 

It seems to be unfeasible task to solve saddle-point equation for the integrand of (fTHj) . 
We therefore take a different approach to the problem, and replace this integrand 
by its asymptotics at ~ such that 

i?j ~ 0{1) , mi + 1712 — is arbitrary , mi — m2 — is fixed . (80) 

Here Ri is classical angular momentum defined in Eq.(|94p. The last condition states 
that the main contribution to the integral comes from diagonal components of SU (2) 
matrix elements, non-diagonal contributions can be treated perturbatively. To be 
precise, we adjust the following three approximations: 

1. We use the following asymptotics for the modified Bessel function 

Ux) = exp(-^n^) (1 + 0{x-')) (81) 

when X ^ 00 and such that n^/x ~ 0(1). It leads to 

e^'^ 1 

l2s{2/3cosiu) = exp(--s^ -/^sin^cj) (l + C>(sin^ cj)) . (82) 

Due to the property (f7^ the remainder is bounded like 0{j3~^). 



2. The integration region over uj is extended to infinity after proper cliange of 
variables (tliis introduces only exponentially small corrections which can be 
properly bounded). 

3. The main approximation concerns the behaviour of matrix elements d^'^iiS) 
in the vicinity 0. The standard results on the asymptotic expansion of 
d^'^iiS) assert that when a; — *• 0, r — >■ cxd the SU{2) matrix elements approach 
the IS0{2) matrix elements pHITij 



where Jm.{x) is the Bessel function. We have found, however that this approxi- 
mation is oversimplified. In particular, it does not lead to the correct Gaussian 
distribution of the dual Boltzmann weight that is clearly seen from the effective 
theory obtained above. Moreover, Eq. ljH^ states that, e.g. the asymptotics of 
all diagonal elements is the same. To get more precise asymptotics one can use 
representation of (i-function in terms of the Jacobi function. The latter can 
be treated with help of the Tricomi expansion into degenerate hypergeometric 
functions as described in ^^1; section 10.14. Finally, we use uniform Tailor 
expansion for the degenerate hypergeometric functions JX^, section 6.13. This 
leads to the following asymptotics 



where all notations are given in ()120|) . It can be easily checked with help of 
Mathematica that the last equation gives better approximation to ^-function 
compared to Eq. ljH^ . especially at large values of m + n and is valid, in fact for 
all values of r. Moreover, it gives very reasonable approximation in rather wide 
region of u, except region ~ 7r/2. In the appendix B we give independent and 
simpler proof of this asymptotic expansion calculating also the first correction 
0{siia^uj) to Eq.()84j). The final result is given in Eq. ()119|) . We term this 
asymptotics as /S'0(2)-like approximation of SU{2) matrix elements. The 
only (but essential) difference from the exact IS0{2) matrix element is the 
presence of sin 6' in the argument of the Bessel function. 

Combining Eqs. (j82p . (j84j) we arrive finally at the following asymptotic represen- 
tation for the one-link integral 



(83) 



(84) 




■mi—m2 



exp I ——5 




(85) 



where we have used notations 




(86) 



a = mi + 1712 ~ ti —t2 



(87) 



Making change of variables sincu = y/^/2(3 in the last integral and extending the 
integration region over y to infinity the last integral can be written as 

^ Jo ^ ^ ^ Jmi-m, \^—^==-y 1 Jt2-ti \^—j==-yj , (88) 



where R and sin 9 are defined in Eq. ()120p . Calculating the last integral we find 

1 



So(/) X C{/3) 5r,r,r exp 



-^{Rl + Rj- 2R1R2 cos 6 1 cos ^2) 

f Rism9iR2sm62\ 
x4( J' (89) 

where k = rui — m2 = ^2 — ^i- Again, we have compared the present asymptotics 
with the exact expression given by Eq. ljHOj) using Matemathica. Suprisingly, if (3 is 
sufficiently large the asymptotic formula gives very reasonable approximation for all 
configurations we have studied. It is interesting to stress that the same asymptotics 
can be obtained directly from ()3U|). where the CG coefficients are replaced by their 
semiclassical limit by Ponzano-Regge ^H] ■ The method presented above has however 
two advantages. First of all, we have computed corrections to the leading behaviour 
(jHH), therefore (9(/3~^)-corrections to the asymptotic expansion of one-link integral 
can be computed in a straightforward manner. Corrections to the Ponzano-Regge 
formula are in general unknown. Secondly, we find that our method is much simpler. 



4 Discussion 

In this article we proposed dual formulation of non-abelian spin models. Our ap- 
proach to the dual transformations is summarized in the beginning of the section 
2.2. Main formulae of the present paper, Eqs. (j^ - (j^ . (pSj) - (j3^ give two ver- 
sions of the dual formulation. To demonstrate the usefulness of our formulations 
we have computed the low-temperature asymptotics of the dual Boltzmann weights 
for two-dimensional model. We reckon that the results for the dual weights provide 
significant simplification of the partition function at large values of (3. Nevertheless, 
even in this case the model is still too complicated to be solved exactly. In our 
next paper [T7] we consider some physical applications of the dual representation. 
In particular, we attempt to give a saddle-point solution of the dual model and 
derive effective ID model for the two-point correlation function. Also, we calculate 
continuum limit of the dual representation and solve classical continuum equations 
of motion. 

Some simple consequences which worth noting can be obtained from ()74p. Let 
us perform a shift of variables in the integrand of ()74|) 

Rk{x) ^ i?,(x)--^Vjfc(/)A(x) . (90) 
27P ; 

This shift allows to extract the main contribution from the sources to the free energy 
which emerges from the Gaussian term in the effective action, i.e. Sn eiven bv Ea.dSSI) 



= C(/?)exp 

where Z{j,f3) is a remaining part of the PF. Here Di{x) and Gur are hnk Green 
functions introduced in In |_17J we show that the following sources 

jkil) = -2i^/j(JTT)hk , leC^y , hl = l 

and jk{l) = 0, otherwise, can be related to the correlation function in the represen- 
tation j. Substituting these sources into one finds after simple analysis 

Tj{R) = Z.(i,/3)/Z.(0,/3) = exp 

where D[R) = 0{lnR), for i? — > cxd. By no means we claim that this result 
establishes power-like decay of the correlation function: in fact, we do not know the 
i?-dependence of the remainder. Therefore, at the moment we can only claim that 
the shift ()9U|) allows to extract the main perturbative contribution. We would like to 
stress that the shift ()90|) is precisely the same as the one used in . Unfortunately, 
we do not know if the methods of 2J can be applied to the partition function (j9ip to 
bound the remainder: the essential difference from the XY model is that even in the 
absence of vortices both full and effective dual actions are complex. Nevertheless, 
we think that such possibility deserves further investigations. 



27/5 



(91) 



j(j + l)DiR) + Oip-'/') 

ll3 



(92) 



A Asymptotics of ln<l>(f; V) 

In this appendix we calculate the asymptotic expansion of the function ln$(^; V^), 
where $(^; V) = $(r, t,u;) is given by Eq. pn)) 

$(r,t,u;) = Y.i~^r P,(t) xU^) • (93) 

A=0 ^ 

Here, Px{t) is the Legendre polynomial and t is given in Eq.()41|). SU{2) matrix is 
taken in the parameterization (|34j) . (|35p. Also, we define and compute the semi- 
classical limit of $(r, t, uj). To this end we introduce classical angular momentum as 



R = 2r + 1 . (94) 

We are interested in the expansion of $ (^; V) when u ^ uniformly valid in all 
variables ^. Because ip and a are compact variables the problem of uniformity 
concerns mainly the value of the variable r. We face two situations: 1) r is fixed 
and 2) r — * 00. In the second case we look for the asymptotic expansion such that 
Ru ~ 0(1). 

As the first step we express the generalized characters entering $(^; V) through 
the associated Legendre functions of the first kind taken on the cut x G [—1,1] 



(2r + l)(2r + l + A)! 

(9.r - W sin H 



(95) 



When r is fixed, we construct Tailor expansion of $(^; V^). It is more convenient to 
expand in powers of y = sin ^ rather than in powers of u. From the representation 
it follows that 



R{R + X)\ 
2r(A + 3/2) [{R-X- 1)!_ 



1/2 



2A + 3 2 ^ 



(96) 

Substituting last expression into Eq. (j9H|l . it is straightforward to find the following 
expansion 

In <l>(r, t, u) = -zyF,{R, t) + ly^F^^R, t) + '-y^F^{R, t) + 0{y^) , (97) 

o o 

where the coefficients of the expansion Fi{R,t) are 



Fi{R,t) = WR^ - 1 



(98) 



F2iR,t) = (3^2-1) 



R' 



1 - - - 4 



(99) 



+ ^{R^-4)-R^ + l- \^/W^VW^ + Vi^^v/^R^KIOO) 

When ^ 1 we expand coefficients Fi in powers of R. This defines the semiclassical 
expansion of $(^; V^) when — > 0, i? — > oo such that -Rcu ~ 0{1) 



(101) 



It would seem that the bound 0{y'^) in ^7\\ and, therefore in ()101|) does not 
hold uniformly in R. It turns out however that the bound is correct both for fixed 
r and for r — oo. To prove this we need an expansion for the associated Legendre 
functions at large r uniformly valid in the neighbourhood of the point u = 0. Such 
an expansion is given by MacDonald's formula ^S], section 3.5 



P;''(coso;) 



/ 1\ ^ 
z/ H — I cos — 

^ 2' 2 



(102) 



X 



J^(a) 



2 ^ 

sin — 
2 



a 1 
-J^+3(a) - J^+2(a) + —J^,+i{a) 



where a = (2z/ + 1) sin ^. The ratio of factorials in (jH^I) is expanded as 



(R + xy. 

_{R-X-1)\_ 



1/2 



^A+1/2 



A(A^ + |A 



0{R- 



(103) 



With these results, and using recursion relations for the Bessel functions to reduce 
order to A + | we obtain the representation for the generalized character of the form 



/f 


R 




9 sin 



1/2 



B{R,X;uj,dh) \h=i Jx+i{cth) . 



(104) 



The function B{R, A; u) has the following asymptotic expansion 

+ 0{sm'^) , (105) 

where a = 2i?sin^. Substituting last expressions into Eq. (j93p we can easily calcu- 
late all the sums and derivatives. This leads to 

ln$(r,t,a;) = -2iRt sin - + —tsin- (106) 

+ sin^ - 3t^ - 1 Rt^ sin - + sin - + C sin^ - ) . 

4V 3 4 4yV4/ 

Re-expanding last formula in powers of sin ^ one sees that Eq. ()106p coincides with 
Eq.JMI). Due to the bound 0(u;^) in dTTl^ . the same bound holds in (ITn?)|) and, 
consequently in ()101|) and in (jOTjl . We thus conclude that the Tailor expansion (jOTjl 
does provide expansion at small lo which is valid both for fixed values of r and for 
r — »• cxo. 



B{R,X;uj,dh 



sin'^ 

1 + 



7 n 9 

4-b^^ + ^)-2h^ 



B Asymptotics of d^^i^uj) 



Here we compute the asymptotic expansion of the SU{2) matrix elements dl^"'{uj) 
in the classical region R = (2r + 1) ^ 1 uniformly valid in the vicinity of the point 
u = for all allowed values of m and n, m — n is held fixed. 

To get such an asymptotics we first present d-function in terms of hypergeometric 
function F = 2-F1 



{s + k + p)\{s + k)\ 



sm 



UJ 



cos ■ 



xF{s + k + 1,—s — p;k + 1; sin — 



where ^„ 



1 if n > m, ^„ 



-1 otherwise, and 



k 



\m — n\, p = \m + n\ 



-ik + p) 



(107) 



(108) 



As is seen from the arguments of the hypergeometric function the infinite series in 
F terminates so that right-hand side of ()107p is polynomial in sin^ ^ 

dTiuj) = Un \A,{x)A,{y)]^ (sin ^ 



r+i(p-fc) 

X E (-1) 



UJ 

COS — 

2 



(sin^ 



r(A; + l + /)/ 



(109) 



where we have introduced the following notations 

r(x 2^ 2) 



Ak(x) 



r(rr-L. Ih-L. 1\ 



l^ ' 



(110) 



x = r + ^{l+p), y = r + ^{l-p) . (112) 

The second step consists in expanding the ratio of Gamma functions. This can be 
done with help of the following formula 

Ux + a) ^a-../'g(zir5(-W)(,)(5-a). + (9(x-^)) , (113) 



T(x + b) \ ^ s\x 

where B'f'^ [x) are the generalized Bernoulli polynomials. Important point concerns 
the large expansion parameter we use. We take not merely classical angular mo- 
mentum r + i but rather quantities x and y defined above. Such a choice gives 
more accurate asymptotics valid in a wider region of parameters. Then, in the 
case of quantity Ti[x^y) the series in flll3|l terminates because a — b = k + I and 
representation ()113|1 becomes exact. It leads to 

J'iix,y) = ixyy+' J2 , .u/li^ ^ 

xBt'^'\lik + l))i?2+^+^)(i(A; + 1) + /) (114) 

what is essentially the desired expansion at large x and y. It follows from the last 
representation that 

:F,ix, y) = ixyy+'[l + hi + k)i---)- 

1 y x Axy 

+ hl + m + k- -l-k-l){-^ + \) + 0{x-\ y-')] . (115) 

24 x'^ 

For v4fc(x) following the same procedure one finds 



i + ^^k{e-i) + o{x 



Akix) =x ^ 



Substituting last expressions into ()lU9p we get after some algebra 



(116) 



c"M-^..(cos|)-x:(7^(0' 



2 \y X J Axy 48 \x'^ y'^ 



^^{i + m + k-mf-i-k-i) (^-^)} 



Here we have extended summation over / to infinity since this introduces corrections 
of the order 0(00'^'^) or less. Recalling now the series representation for the Bessel 
function 




(118) 



we can easily sum up all series in the last formula. Finally, we arrive at the following 
asymptotic expansion for d-function 

= UnUkit) + ^[Jk{t) + :r^(l - 2cos2 0){Jk-i{t) - 
4 3 sm 6 

Jk-,i{t)) - ^^^ff [{k + l)J,-2{t) -{k- l)J,^2{m + 0(sin^ |)} . (119) 

We have introduced here the following notations: 

i? = 2r + l, cos0^^, 6 = sin2-, t ^ RsmOVb . (120) 
R 2 
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